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Abstract 


We study the thermodynamics and non-relativistic hydrodynamics of the 
holographic fluid on a finite cutoff surface in the Gauss-Bonnet gravity. It 
is shown that the isentropic flow of the fluid is equivalent to a radial com- 
ponent of gravitational field equations. We use the non-relativistic fluid ex- 
pansion method to study the Einstein-Maxwell-dilaton system with a nega- 
tive cosmological constant, and obtain the holographic incompressible forced 
Navier-Stokes equations of the dual fluid at AdS boundary and at a finite cutoff 
surface, respectively. The concrete forms of external forces are given. 
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1 Introduction 


Using the fluid/gravity correspondence [1, 2], it was suggested that the problem of Navier- 
Stokes (NS) turbulence might be mapped to a problem in general relativity [3, 4], with 
different scales appearing in turbulent phenomena corresponding to different radii in the 
dual geometry [5]. Thus, with the holographic Wilsonian renormalization group (RG) 
approach [6, 7, 8], it is also interesting to study the holographic hydrodynamics at a finite 
cutoff surface directly [5, 6, 9]. In this study the timelike cutoff surface plays an important 
role so that the asymptotical behavior of spacetime geometry becomes irrelevant. This 
provides a potential approach to study holography in asymptotically flat spacetimes. 

It was proposed in [10] that there is a mathematically precise relationship between the 
unforced incompressible NS equations in p+ 1 dimensions and vacuum Einstein equations 
in p+2 dimensions and their solutions. The dual geometry has an intrinsically flat timelike 
cutoff surface whose extrinsic curvature is identified as the stress energy tensor of the dual 
fluid. This relationship has been further developed in the literature such as [11, 12, 13, 14, 
15, 16, 17, 18]. In [10], a gravitational shock wave was introduced to stir the fluid and then 
left to evolve according to the unforced NS equations. In order to study the stationary 
NS turbulence, it is better to introduce the external random source fields [19, 20]. For 
example, a dilaton field was added to the bulk gravity in [3]. In this case, a perturbed 
gravity solution with a slowly varying dilaton leads to a slowly varying force term in NS 
equations. Another example is [4], where the force terms come from the slow variation of 
the boundary background in the holographic context. It turns out that a simple forced 
steady state shear solution to the forced NS equations becomes unstable and may translate 
into turbulence at high enough Reynolds number. 

In the fluid/gravity correspondence, the derivative expansion method [1] leads to the 
equations of motion of relativistic fluid, and the equations reduce to the NS equations in 
the non-relativistic limit [4, 21]. Actually, the NS equations can also be obtained via taking 
the non-relativistic expansion method directly [10, 11]. This method can be used not only 
for dual fluid at the AdS boundary, but also at a finite cutoff surface with flat induced 
metric [12]. In the latter case, the bulk geometry is not required to be asymptotically 
AdS. While in the asymptotically AdS case, the perturbed gravity solutions with Dirichlet 
condition at the cutoff surface can be mapped to the perturbed gravity solutions without 
the cutoff surface [22, 23]. To study the holographic fluid with external forces in non- 
relativistic limit, instead of using the induced metric perturbations method proposed in 
[4, 23], we keep the induced metric flat and add external matter fields in the bulk as the 
source terms to the dual fluids. 


It was shown in [5] that the radial component of Einstein equations is equivalent to the 
isentropy equation of the dual fluid. Using a more generic static metric in some sense, we 
show that this statement is also true in the Gauss-Bonnet gravity. Using the non-relativistic 
expansion method, we give the procedure to obtain the perturbed solutions of the bulk 
gravity with a finite cutoff surface up to the second order of non-relativistic expansion 
parameter, and corresponding NS equations of the dual fluid at the cutoff surface. It 
turns out that the shear viscosity over the entropy density of the fluid dual to the Gauss- 
Bonnet gravity does not run with the cutoff surface. This part acts as the service to 
introduce the non-relativistic expansion method. In this paper we mainly focus on an 
Ejinstein-Maxwell-dilaton system with a negative cosmological constant and pay attention 
to the external force terms in the dual non-relativistic fluid. The external force terms come 
from the Maxwell field and dilaton field in the system. Note that the forced fluids at the 
AdS boundary have been discussed in the Einstein-dilaton system [3] and Einstein- Maxwell 
system [24, 25, 26], respectively. The perturbed solutions have been obtained to the second 
order of the derivative expansion. We consider the Einstein-Maxwell-dilaton system and 
obtain the perturbed solutions up to the second order of the non-relativistic expansion 
with/without the cutoff surface. Associated forced NS equations are also derived in both 
cases. The concrete expressions of external forces of the dual fluid, which come from the 
Maxwell field and dilaton field, are given. The results show that the Reynolds number of 
the dual fluid becomes larger and larger when the cutoff surface approaches the horizon of 
the background black branes. 

This paper is organized as follows. In Section 2, we start with a generic static black 
brane metric and obtain the perturbed solutions with a finite cutoff surface in the Gauss- 
Bonnet gravity by using the non-relativistic expansion method. This section acts as to 
fix the notations in this paper and to introduce the non-relativistic expansion method. In 
Section 3, we consider the Einstein-Maxwell-dilaton system with a negative cosmological 
constant and discuss the dynamics of dual fluid on the AdS boundary in non-relativistic 
limit. We generalize the discussions to the case with a finite cutoff surface in Section 4. 
The conclusions are given in Section 5. 


2 Holographic fluid at a finite cutoff surface 


This section is a generalization of discussions in [12] on the thermodynamics and hydro- 
dynamics of dual fluid at a finite cutoff surface. Slight differently, we start with a more 
general static metric and work in the intrinsic coordinates on the cutoff surface directly. 


2.1 Thermodynamics of the dual fluid at the cutoff surface 


To study a fluid in a (p + 1)-dimensional flat spacetime, we consider the generic (p + 2)- 
dimensional static background 


dees = —gu(r)dt? + Grr(r)dr? + Gou(r)disdx'da? , Si Jep = L2 cee, (1) 


where the metric components are functions of radial coordinate r only. We assume the 
metric has a well-defined event horizon at r = rp, where gulr) has a first order zero 
gu(Th) = 0, and g(r) has a first order pole g7 (ra) = 0 [6]. For example, the ingoing- 
Rindler form of flat spacetime and the black p-brane solutions in asymptotically AdS 
spacetime have the form [5]. Using the Eddington-Finkelstein coordinate 7 defined by 


dr = dt + \/grr(r)/gu(r)dr [6], we can rewrite the metric (1) as 
dasa = 24/ gu(r)grr(r)drdr — gulr)dr? + Jza(r)óijdz'dat, (2) 


which has the translational invariance in 7 and x’ directions. We can always introduce 
a finite cutoff surface X, at r = re outside the horizon with the intrinsic coordinates, 
T° ~ (7, 7°), as 


= # = y gulro) T, Psat) T$, {a,b...} = 0, 1,2, ..., p. (3) 


Then the associated bulk metric (2) becomes 


GAY 


das = grr(r)dr? + Yalr) [d7 + V(r) d2dr] az” + N (r)ôżdr] ; (4) 


where Ya(r) and N (r) are given by 


Yar(r)de*de” = — GAT) g2 q Geel) gs gai (pr) = N Iar) Irr (r) (5) 


Julre) Irale) Julr) 
And the induced metric with intrinsic coordinates on X, is simply given as 
d52,1 = Yas(Te)dk"dz” = ad7 d7? = —d7? + dydd. (6) 


The Brown-York stress energy tensor 7,2” evaluated on the cutoff hypersurface £, is 
proposed as the stress energy tensor of the dual fluid [10]. It has a close relation with 
the extrinsic curvature tensor Ka, = sL wYab(7)|r=r, Of Ue, where Lx is the Lie derivative 
along the unit normal N“ of the hypersurface. To study the thermodynamic properties of 
the dual fluid at the cutoff surface, we begin with the re-scaled metric (4) with a Killing 
horizon located at r = ra. The local temperature T, (re) on Xe, which is identified as the 
temperature of dual fluid [5], meets the Tolman relation with Hawking temperature Ty of 
the black brane metric (2) 


= T 1 
Ieo EE = Ty = lim gutr) 


AF] d O ”) 


To discuss the local entropy density of the dual fluid, we consider a quotient of the general 
geometry (2) under shift of 2"[5], 2’ ~ x2’ +€pn', with a characteristic length lo and n’ € Z. 
Equivalently, using metric (4), the spatial RP on X. turns out to be a p-tours T? with 
r--dependent volume 


Volre) = 922(re)Vo, Vo= Gh, ar + Lo Gaal(re)n’. (8) 
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As the total Bekenstein-Hawking entropy S is fixed, we can identify the dual fluid’s entropy 
density as 


ie S _ 1 gee (rn) g — Vol) 
ee Vlr) AGp+2 gh? (re) AGp +2 


(9) 


where S' is the Bekenstein-Hawking entropy of the black brane. 
As a calculation example, we consider the Einstein-Hilbert action with a non-positive 
cosmological constant A < 0, as well as the Gauss-Bonnet term Leg = R? — 4RapRA? + 


RapcpR4°©°, and the contribution from matters Lm, 
1 
[= o) d+? 4,/—9 (R-2A+aLcp + Lm), (10) 
167Gy+2 


where a is the Gauss-Bonnet coefficient with the same dimension as square of length. 
The Gauss-Bonnet term is a topological invariant when p = 2, and will modify Einstein 
gravitational field equations when p > 3 [27, 29]. We will generally demand p > 2 and set 
167G')42 = 1 throughout this paper t. Varying the action (10) with respect to the induced 
metric y~% with appropriate surface terms [31, 32, 33, 34], we can get the Brown-York stress 
energy tensor yaa evaluated on the cutoff hypersurface Xe of the Gauss-Bonnet gravity. 
Because of the flatness of the cutoff surface, its final form is given by [12] 


Te =2 [Kias — Ku — 2a (3.00 ~~ Tia) F Cis] , (11) 
with i 

J y= ZOK Kock + Keg“ hg — 2K K Ku- K?Ko) , (12) 
where K is the trace of extrinsic curvature tensor Ka of Xe, and C is an ambiguous 


constant. Substituting the metric (4) into the Brown-York stress tensor (11), we can get 
the stress tensor of the holographic fluid dual to the Gauss-Bonnet gravity as 


FB" ditdx? = 2 (r.) — 2C| dz? + [ðe (re) + 2¢| S; dž'dž, (13) 


with constant pressure pg (re) + 2C and energy density pare) = 2C’, where 


Pale) = G(r) Pa (re); Üg (re) = Palro) + Bete) = G” (Te) tig (re), 
P Mele) x =A (2 a . (14) 


Grr (Te) Jasle) Grr (Te) a 


Ďe (re) =~ Git(Te) eal Te) 


The sum w,(r-) is independent of the constant Č. The Gauss-Bonnet term corrections, 
compared to Einstein gravity with Ag = (p — 1)(p — 2)a, are given by 


Ste \ pe ÀG Daae) ? Swf) 1] ÀG ato) ? 
MAR ae (5). Wr) = 1 oe mi (5) 


'Most of the results in this paper have been checked by Mathematica up to p = 5. 
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Using the above results, in general we have the thermodynamic relation 


Ires) = We (Te) — 0, (re), Q, (re) = Su (re)MG) (Te), (16) 


where m = 1,2,..., and every gh (re) corresponds to different kinds of charge density with 


corresponding conjugate chemical potential me (re). Note that in the Gauss-Bonnet theory, 
the area formula of black hole entropy no longer holds; but for the black brane solution 
with Ricci flat horizon considered in the present paper, the area formula for the black brane 
entropy still works [28, 29, 49]. Thus the total entropy is given by ° 


S= s4(reVo(re) = > [erd —Qo(re)] as = 0. (17) 


This isentropic equation 0,.,.5 = 0 can be considered as either an adiabatic thermodynamic 
process of the dual fluid or a holographic renormalization group flow [5]. Varying the 
action (10) with respect to metric ĝ4#, we can obtain the equations of motion for the 
Gauss-Bonnet gravity in the re-scaled (p+2)-dimensional bulk coordinates 


Wap(r) = Eapg(r) + aH ap(r) — Tag(r) =), {A, B...}= {r,7,z'} (18) 


with the definitions 


~ 5 Tan 7 z ôL z 
Ean(r) = Rag — ze gas + AJAB, Tag(r) = 2 an + JABLM, 
k (19) 


z ` p 7 ` = = D L 
Hap(r) = (RacporRg  — 2RaceoR” — 2RacR° , + RRap) — 3JABLGB, 


where T4g(r) is the stress energy tensor of bulk matters. If we assume the metric (4) solves 
Eqs. (18), the nonzero components of the stress energy tensor are T(r), T(r), Tza(r) = 
Txizi(r), and T;,(r) = T.z(r) = N(r)T+(r), respectively. Following the procedure in [5] and 
using metric (4), we can introduce an arbitrary null vector ¢ tangent to the hypersurface 
r = rą with time component 0;, such as 


C404=0;-On Gane“C?=0, C Na =n, (20) 
with N the unit normal of the cutoff surface. It could be checked that 


a baer alr) p D Irr (re) gu (re) 2 CIC” |Ean(re) +aHap(re)] . (21) 


When the &, (re) term satisfy the following relation 3 


Vo(Te) e _ Volre) = ro CACB Ty alr 
ar, 7 “26, (r9| = HO Sirid) PFa. (22) 


In the fluid/gravity correspondence, the corresponding Wald entropy current of curvature-squared 
gravity theories is studied in [38]. 

3For an example see Section 4 of this paper or [5]. For the vacuum case without the matter source in 
(10), this relation always holds. 


Using (17),(18) and the above two equations, we have 


0,8 = BE Sg Credaulre) 2 ACP Wale) (23) 


AS \/Grr(Te)gte(Te), Vp(re) and Ty have neither zero nor pole even when re = rp, we arrive 
at 


ðS =0 — > C40?Waa(r.) = 0, (24) 


which implies the equivalence between the isentropy of the RG flow and a radial gravita- 
tional field equation of the Gauss-Bonnet gravity. Thus we have generalized the result in 
[5] to the case of Gauss-Bonnet gravity. 


2.2 Incompressible Navier-Stokes equations from gravity 


Keeping the intrinsic metric (6) of X, flat, we can take two linear diffeomorphisms based 
on metric (4): a transformation of the radial r and a Lorentz boost in (7, £') coordinates. 
The first is the transformation of r and the associated re-scaling of (7, £*) as 


Jtt (re) ği at Ira (72) 
gulk(r) SN deke) 


where k(r) is a linear function of r, whose form would be chosen via the global symmetry 
of geometry (4). The metric (4) is then transformed into 


rok(r), TOT (25) 


nO pgr HO gp | Geol?) 


Jtt (Fe) Jtt (Fe) Ira (Fe) 


dsp, > dêz = 2 idx dat. (26) 
where we have introduced the re-scaled coordinate f = k(r) and the notation fe = k(re). 
We will work in the (r, £*) coordinates directly in this paper. 

The other diffeomorphism is the Lorentz boost with a constant boost parameter 8; in 
the 7° = (7, Z') coordinates 


Fig, BALE, A= (~s dt (y= 1) (27) 


where we also have defined the (p + 1)-velocity 
=i ee 
Ue = 7 el bi), Ve = (1 = 8°) i ’ B? = pib’ = dib P. (28) 
As these parameters are all constants, the associated transformations can be expressed as 
d? + —ti,dt*, dž > M dł, ðpdt'd > 6, ni nidzdz’ = Pydz°dz’, (29) 
with the projection operator P= Aap + tatty = inini. Then the metric (26) becomes 


d8? 3 > d3? o = GF) dF? + Fan(*) jaz" + N (ryiitar [az +N(*Aabar], (30) 


T 


where we have used the ADM-like decomposition at the constant 7 surface, 


6 A  Iulf) anu Gal?) > Y F= = Gt (Te) Grr (F) 
ae gegyree An Pes 


After the two diffeomorphisms, the metric (4) and (30) still solve the same gravity field 
equations. The cutoff surface r = re is equivalent to 7 = fe, we will firstly work with the 
new radial variable 7 , and then transform back to r via r = k7'(7). 

In general, the Brown-York stress energy tensor on the cutoff surface r = re corre- 
sponding to metric (30) will turn out to take the form 


Ta Oe E P(Tc)tatly + Dre) Pav, w(Te) = Ge) E5 Ge (32) 


which could be identified as the stress tensor of an ideal relativistic fluid Lo with 
(p + 1)-velocity ŭa in flat space-time. In the Gauss-Bonnet gravity case (11), the cutoff- 
dependent energy density P(r.) = pg (Fe) —2C and pressure f(r.) = Dale) +2C. In general 
here C is an unfixed constant, but to obtain a finite result when the cutoff surface goes to 
the AdS boundary, Ĉ can be fixed [12]. 

Let us pause to recall the hydrodynamical description of microscopic field dynamics 
in flat spacetime. It applies when the correlation length of the fluid leor is much smaller 
than the characteristic scale L of variations of the macroscopic fields [35]. Via dimensional 
analysis, leor ~ 1/T, and 1/L ~ Oza, where T, is the characteristic temperature of the fluid 
and za would act as the coordinate-dependent parameters, such as T(%) and t7(%) *. One 
introduces the dimensionless Knudsen number Kn = leor/ L ~ 43, ~€ <1 to expand the 
stress energy tensor of relativistic fluid in flat background ‘ap, 


Tav(®) = Ty (2), Ty (#) = WB) tia) tio £) + B(@) fran, (33) 


and TAY (z) ~ (Kn)”. To take the non-relativistic limit as in [35], we recover the speed of 
light c in (33) and introduce the normalized pressure P(%) as 


fia(@) =, (8) (-1,B(B)/c), 7,(@) = (1- B(@)/2)7?, BE) = BEE), a 
P J ï) TE) l ÕPE alz) O,P(z) (34 
rað = rE ar @) a, —— ~ — ~ — en i 


When c — oo, the energy-momentum conservation equations of the ideal fluid, 627 (2) = 
0, leads to the non-relativistic incompressible Euler’s equations, 


O,P(#) + 4,8;(2) + B7(#)0;8(2) =0, p (t) =0. (35) 


“We have used (#) to denote the function arguments (°), and would use 0, to denote Aja. 


Instead of the c — oo limit, if we assume the small velocity parameter to take the same 
limit as the Knudsen number, ie. 3/c ~ Kn ~ e€, it would be equivalent to set c = 1 in 
(34) with the following scalings 


ðe Orne, Bl#)~e Pee. (36) 


These non-relativistic scalings are named as the BMW limit [4, 23]. To obtain the forced 
NS equations, the dissipative part of the stress energy tensor (33) is required 


Tw(@) =TS O+T'@, TE @ =THO+TY@+... 87) 


a 


For example, in the Landau frame é*(#)7,2"**(#) = 0, the first order dissipative components 
could be written as 


Top’ (E) = ~2ii(@)Bao(@) — CO E)PavlE), Posl) = an + tial Bil), 
P 


: 7 7 (38) 
Öl) = Pa (£)P 3 (2) âmin (7) — 


where #(%) is the kinetic shear viscosity and (č) is the bulk viscosity, the latter vanishes for 
conformal fluids. Usually, they behave the same as local temperature T(Z) in (34). The first 
order dissipative hydrodynamics satisfies the dynamical equations 0% as (z) + TP (E) = 
AS) if some external source term appears. In the non-relativistic limit (36), if we further 
assume f;(%) ~ e, fz(@) ~ ef, it would lead to the non-relativistic forced incompressible 
NS equations at order e, 


ÕP (t) + 8,8;(@) + Bi (EÂ LE + ael = f"(®, AlE =o, (89) 


where (%) = n(Z)/w(%) is the dynamical shear viscosity and fla) = f,(@)/w(%). The 
NS equations have the scaling symmetry (36), as shown in [4, 10]. 

Next we will derive the incompressible NS equations from the gravity side. To get 
the dissipative part of the dual fluids, we need to perturb the geometry (30) by using 
either the linear response method (eg. [6]) or the perturbative expansion method developed 
in [1], where a perturbative procedure to solve Einstein’s equations order by order in the 
boundary derivative expansion was proposed. In the case with a finite cutoff surface, 
we also can employ the procedure under the non-relativistic limit [11, 12]. Regarding 
the transformation parameters in linear function 7 = k(r) and the p-velocity 3° in the 
boost transformation as functions of the hypersurface coordinates (7,Z*), we can solve 
the gravitational equations order by order in the non-relativistic perturbative expansion. 
Based on (36), these transformation parameters are also regarded as small quantities with 
appropriate scaling symmetry as 


ð, ~ el, O~6(#) +e, ð ~ dk(r,%) ~ P(z) ~ e, (40) 


where 6k(r) = ř—r ~ e? would provide the pressure perturbation. Using the formal Taylor 
expansion such as gulf) = gu(r)+g1,(r)Ok(r), and define gz,(r) = V/git(T)Grr(T) /gut(Te), we 
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can consider the metric (30) with coordinate-dependent parameters as d5)» and expand 


it up to order e, ° 


d5(o) = + 2Ger(r)dFdr + |- gulr) gr Sel) 5 ari aw 
(0) =+ 25r) wod rd" 
( 


= Mr) Baa + 2 [IE — See] awar 


4 Geer) 82(B azar — | LOO) _ Ilr) | (5205) a2? + 8.(2)8, (Baz aR 
Ger (r)B?(#)dFd — — (6°(a)d7* + B,(#)B;(&)dz'dz’) 
jee Irv (T) ay I(T) Jal re) 5 = gulr) = 
rl E j ao omg )— Jul a kro) (a rr nae n AO ‘ 
deie r — Izel) Joa(T)OK(r) _ Gale) r Z, gi 
aa E Eo a a J Ga 
+ O(e), (41) 


where 6k(r) and 3;(%) and are all functions of Z° now, we will omit this notation (%) hence- 
forth. The coordinate-dependent metric (41) is no longer a diffeomorphism of metric (30). 
Under the perturbative expansion; it turns out that (41) only solves we same equations 
of motion of gravity up to et. To solve the equations of motion up to €?, some constraint 
equations and new e terms to the bulk metric are needed at ihi order. With the 
same way, the equations of motion of the system can be solved order by order in the non- 
relativistic expansion parameter e [11]. And the corresponding Brown-York stress energy 
tensor which is identified as the dual fluid’s stress energy tenor, can also be obtained at 
the desired order. In this paper we solve the equations of motion up to order e°. It turns 
out that the non-dissipative part is still given by (32) in the non-relativistic limit. Up to 
e€, we have 


Ty didi? = [A(r) + Ü, (re) B?] dF? — 2tdo(r.) bid d7 
ab 0 


` pok (42) 
+ [B(r'e) diz + (re) Ab] dz dz + Olè), 


where Ñ, (re) = A (re) + Bo (re), and 
Ple) = Bo(Fe) S Palro) + Pi(e)Ok(re), Pre) =Do(Fe) = Po lre) + Bi (Te)OK(re). (43) 


The equations of motion of the bulk matters should be solved in the same procedure. Once 
we get the solutions of both gravity and matters up to order e?, the constraint equations of 
gravity at order e are just the forced incompressible NS equations (39) of the dual fluid. 
As a calculation example to perturb the geometry (30), we again assume the metric 
solves the equations (18) of Gauss-Bonnet gravity with a non-positive cosmology constant 


>In this paper, we use the scripts (0), (1), ... to denote the order in the derivative expansion, and scripts 
(€), (€)... to denote the order of the non-relativistic expansion parameter <€. 
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A. The metric (41) leads to we, = 0, and the constraint equation at order e° turns out 
to be A 7 

INOW) = -ETO = a, (r-)d,8' = 0, (44) 
where wÙ, (re) = W,(r-) is nonzero outside the horizon. Thus it leads to the incompressible 
condition 0;3' = 0 of the dual fluid at the cutoff surface. 

To solve gravitational field equations at order e°, we need to add corrections to the 
metric (41). It was shown in [1] that due to the spatial SO(p) rotation symmetry of the 
black brane background, one has the decoupled equations of SO(p) scalar, vector and 
traceless tensor perturbations. For example, if we turn off the tensor perturbations of the 
matter stress energy tensor Tag(r), the tensor corrections to the metric, at order €°, are 


F- TT È ~ A 1 A 
d3) = oF (nr) Sel Suelo) 5 Dij dz' di) + a Tij = dubi) = ~ 55,0, 8". (45) 
Gux(Te) p 
Then the traceless tensor parts of the gravitational field equations WwW) = 0 lead to a 


ij 
second-order ordinary differential equation 


gal) m Vary) [, _ Golrn) 
(ee io) en 0 Fi) = fo Tea) "ace ae ay, oo 


where two integration constants have been fixed and 


joutr) g(r) okel) 


—1)/2 
(p—1)/ (r) 


G(r) = 4s l—alp—2 
Peed eee 


(47) 


In (46), one integration constant is the upper limit of the integration, which is chosen to 
be re via the Dirichlet boundary condition at the cutoff surface. The other integration 
constant is chosen to be G,(r;,) to cancel out the first order zero in the denominator 


VJ gut(Tr)// Grr(Tn) = 0 of the integrand in (46), where 


i= lim — „p/2 rh = ada A. = A(p — 2)nTy IoalTh) 
Oln) = Jim G0) = BP- aha), Aa = ETE fat 


; (48) 


and 4/gulrn)grr(rh) is a finite constant due to our definition. The Hawking temperature 
Ty is given in (7), so that A, is a constant determined by the metric (1) at the horizon. 
With the perturbed metric d?o y+ dst), corresponding Brown-York tensor becomes 


qo) + J , where the symmetric ees components of 7) are 


JAD g Golts) | g VITO gry] g,, — 2 Galt) 5. = anle): 49 
7 dE ‘ll ee “eee ee 


Here 7(r,) is defined as the kinetic shear viscosity of fluid dual to the Gauss-Bonnet gravity 
with the geometry (1). Using the entropy density s,(r.) given in (9) of the fluid at the 
cutoff surface, we have 


“iy — Salts), lr) in, 
re) gee (re) 7 So (rc) T i Aal (50) 


One can see that the shear viscosity over entropy density of dual fluid does not run with 
the cutoff surface. This match with the results in [6] and [36], where a hypersurface in the 
given solution is introduced. Here we impose the Dirichlet boundary conditions [5, 12] and 
find the corresponding perturbed bulk solutions with a flat induced hypersurface. Using 
the formula (50), one can easily obtain the ratio of shear viscosity over entropy density 
for some fluid once the dual gravity background is known. Take an example, using the 
charged AdS Gauss-Bonnet black brane background solutions, one has [15, 30] 


2 2 p+1 2p 
Or 4 aT Arg Th Th 
Cal?) a g2’ gulr) — Irr (r) — ye 1 = 1 = 2 h E (1 + Q) + | , 
Are) 1 Ag (pt+1 2 
S E a E = (p— 1)\(p—- 2 1 
w) ih- (2e) de=@- 10-2), G1) 


where a is the Gauss-Bonnet coefficient, and £ is the radius of AdS spacetime. When 
L — œ, the negative cosmological constant A = —p(p + 1)/2¢? — 0, the black brane 
solution will degenerate into the ingoing Rindler space via some coordinate redefinition [5], 
and the Gauss-Bonnet term will not contribute to the shear viscosity any more. This has 
been shown via using the ingoing Rindler metric of flat spacetime [37]. Here our formula 
in (48) leads to this result directly because grz(r) = 1 > Ag =0. 


3 Holographic forced fluid at AdS boundary 


It is interesting to give the concrete expressions of some external forces of holographic fluid. 
For this, let us consider the Einstein-Maxwell action with a dilaton field © and appropriate 
boundary terms and counter terms as follows. 


1 1 1 
= ——— | d’**z,/—g | R- 2A — —F? — ~(V®)? | + Sha + Set, 52 
167Gy 42 J INH ( g” gV) + Soa + Ses. (52) 
where the negative cosmological constant A = —p(p + 1)/2¢?, the Maxwell field Fun = 
2V mAwn), and F? = Fun F”, (VE)? = Vm®V S , where {M, N,...} run over the bulk 
coordinates. In what follows, 167G,,2 = 1 will be adopted. From the action we can get 
the gravitational field equations 
1 


; 1 
Run — NR + Agmn = TMN + Tun: (53) 
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where TG), = FypFy? — 2° F?, TẸ), = Vu ®V yO — 24 (V6)?, and from which one has 


#1 1 1 1 
Wun = Run + Py Imn = gm POn® -73 (Puer = gouv?) = 0. (54) 
The equations of motion of the Maxwell field and the dilaton field are 
Wy =VuFX=0, W, =V?O=0. (55) 


The system has a class of exact solutions with p + 2 parameters as 


2 
ds? = —2u,da"dr + a [V(r)uju,detde” + Pry] detdx”, 


A= -U(r)u,dx", 
= ĝo. 


(56) 


where {j1, v, ...} run over the boundary coordinates. This is a charged boosted black brane 
solution in Eddington-Finkelstein coordinates with a constant dilaton o, and a constant 
unit normalized velocity u,, 


ai o 
Up = W%w(—1, vi), Ww = (1 — v’) E vu? = yv' = b,j0'v". (57) 
and 7“”u,u, = —1. In the metric of the solutions (56), 
me Êe oi ar? 
Viper ap = a Oe (58) 
Here the notations in [25] have been used, 
me l 
V(r) =1-M+Q? =0, M= 1? ==. (59) 
ri Th 


where r = rp is the horizon location of the charged black brane solution. The Hawking 
temperature of the black brane is given by 


1 fr2Vv(r)] ay Th 
Ty = — |—— = 1)M -2 — 
H ÅT | (2 | an [(p + ) p Q ] Arl? (60) 
In the gauge field part, we have chosen the gauge that A, = 0 and 
p-1 
_ Th — q 2p 
U(r) = po (% = t) ; Ho = yet pl’ (61) 


where bouo is the value of the gauge potential at the AdS boundary. Usually bọ = 1 is 
required to get a well defined gauge field at the horizon [39, 40], but this is not necessary 
in the perturbative fluid/gravity correspondence [24]. 
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3.1 Charged fluid in non-relativistic limit 


To perturb the metric (56), we take the same procedure via regarding the parameters 
in the metric and gauge field as functions of boundary coordinates u, —> up lx), Th > 
ra(1 + P(x)) ê. Without loss of generality, we take the associated replacements as m — 
m(1+P(a))?** and q > q(1+P(z))” to keep M and Q as two constants. The scalar 
field is replaced as ® + ¢(a"). Then the solutions (56) will not solve the equations of 
motion (55) and (54) any more. Using the so-called BMW limit at the boundary [4, 10] 


Oy ~ e, ĝi ~ vi ~ Ah ~ et, 3x Pre, (62) 


where we have added an assumption that jọ ~ et, we can solve the equations of motion 
order by order in the small parameter e [11, 12]. Actually the three parameters (m, q, ra) 
relate to each other via V(r;,) = 0, so that two of them are independent parameters, for 
example, one can take q and r, as two independent parameters. We focus on the forced NS 
equations in this paper, so 6m ~ ôq ~ orp, ~ e have been assumed in the non-relativistic 
limit. In the boundary derivative expansion of the relativistic fluid, the solutions of the 
Einstein-dilaton system [3], and the Einstein-Maxwell system [24, 25] have been obtained 
up to second order of the derivative expansion parameter. We can extract the relevant 
terms up to the second order in our non-relativistic expansion parameter e. To compare 
with their results, we write our the non-relativistic perturbative solutions up to e° in a 
covariation form, 


2 


ds? = — 2u,,(x)da"dr + 5 |-V (F) tlt tle) Pela onan” 
2 
r V 
+ B [PKR (r) Ogu, + 20F (row + OF s(r) (06) |] dz” dz”, (63) 
A =A,(r, x)dz” = [AW (r, £) + Al) (x)] dt"; 


® =d(x) + LF (r) (u3 p) + C Falr) (0 3 o). 
Although higher order terms such as u,u; ~ e° are reserved, they would not contribute to 
the calculation up to order €°. In (63), Pus (£) = Qu + Up(x)u,(x), F = [1 — P(x)] r, and 
A\ (r, £) = — [1 + P(x)] U(F)u, is the perturbed and boosted potential. Ale?) is an extra 
electromagnetic field, which would provide extra forces and FKP = Jopa are chosen to 


meet with the appropriate non-relativistic expansion in [5], Pe) ~ e, ne ~ e°, as well 


as the velocity u, = (—1—v?/2, vi). This solution solves the equations of motion (55) and 
(54) up to order e°, if the following expressions are provided 


a (y? = Th) Piw 
=f 2 d y = uun = 0, 64 
Po J yeV(y OMe CH) y ee) 


°We briefly denote the function variables (x) as (x) and some of them in the following equations will 
be ignored. 
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where 6 = n” 0 puy = P'”0,,uy, for the dilaton shear perturbations 


œ fa =i 
roA S C oo = Aub Po 6) 
and for the dilaton scalar perturbations 
(2 
Ke) = soa 04 = PO G06. (66) 


We can generalize the dual boundary fluid in [3, 25] to the (p+1)-dimensional case with 
p > 2, where the stress tensor 7,” and current J” are given by 


p+1 9 9 a 7 
TH = iim |2(Kagh®®ht — KH) — Pye 4 2 Gu —_§ (ayy , 
roo (PT L p= 1 p= 1 
prt 1 1 (67) 
H} ur ^p pe = Pp ~ y DIr —  (N\2he 
TJT lim pri Nee Ge =R 5 ith, (Gal, = D! D,® zP) hfs 


where N is the outward pointing unit normal of the regulated boundary, D is the covariant 
derivative associated with the boundary metric h#, and the corresponding Einstein tensor 
G#, as well as the effective counter tensor (Gs) of the dilaton field [41, 42, 43]. Substituting 
our solutions (63) in (67), the Brown-York stress tensor and the induced current at the 
AdS boundary are given by 


i ae ie 1 7 

V Ha, Y Y] 9 Lh pW ; pw 

TH = Mal + (P+ DPI[@ +1) wu? +n] — 27 ot — TR (oo) m 
PE i r? qy/2p(p — 1) 

Th = GEM Sng, ng =(L+pP)no, w=gpU t) = 


where nq) is the induced charge density. Note that here (a4)"” is included in the stress 
tensor (68). These terms relating to the dilaton field were considered as the dissipative 
parts of the dual fluids’ stress tensor in [3], or the perturbations in the membrane paradigm 
[6]. In this paper, however, we focus on the forced NS equations and move them to the 
right-hand side of the constraint equations as additional sources, similar to what has been 
done for the boundary metric perturbations [4, 23]. Then we can redefine the stress tensor 
of the dual fluids with first order dissipative term 


yi 
w Nuw Th 
T = Wo [1 + (p+ 1)P] py + 2e — 210 Cw, Wo = (p+ 1) M a (69) 


where 7 is the holographic shear viscosity of the dual fluid. With the corresponding 
entropy density Sọ, we have 


to = > so = w oe p (70) 


and the thermodynamic relation T7759 = Wo—NofMo Still holds. Furthermore, we can redefine 
the following stress tensor by moving out the terms associated with the boundary chemical 
potential 


T9 = Tyso |1 + (p + 1) P] ur + — — 2 M Cuv, (71) 


This form will be used later. 


3.2 Forced Navier-Stokes equations 


Using Gauss-Codazzi-Mainardi relations near the boundary [3], we can obtain the con- 
straint equations of gravity 


2D4 (Kenh Ph — Kg) = —2RophGN” = -DBV pË N? — FopF” phg N?, (72) 


where hag, N^ and Kap are respectively the induced metric, the outward pointing unit 
normal and the extrinsic curvature of the constant r hypersurface, with the covariant 
derivative D. When r — oo, transforming these results into the boundary metric nuy of 
the dual fluid, we have 


er 


0,7? = — lim 


roo P+1 


[(v°vee + Do) D,® + N°FeoF,?| . (73) 


Introduce the gauge field tensor F, ie = 23 AG” which comes from the boundary chemical 


potential Ave) = —bouo(1+P)u, , and denote the background gauge field tenor as F; (ba) = 
20, AGI, where A?” = ACO + A) the (p +1) constraint equations up to eè become 
ir 1 a 


OT w = —OvP $ (u" One) + pole 


. (0°0,8)] - THEY. (74) 


where FË” = FS”) + FÈ”, and the constraint equation of Maxwell field reads ô" J, = 
N(q)O"u, = 0. Using the redefined stress energy tensor (69), we get the forced NS equations 
at order e, 


OMT = fH 4 fO, fo = -ngut Fes) (75) 


Ww > 


where f9 is just the Lorentz force of the charged fluid and f is the external force from 
the dilaton field up to e 


AL = — xo (W8 00") + & ("08,00 A gs2 An 
= — U — =. = — ——_.. 
The Lorentz force up to eê can also be divided into two parts as 
BP SLO + 1, PO = gut Fn”, fe) = maul). (77) 
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In addition, we notice that the following relation holds at order e° 


BHT) — TD)]bp = FOO A orp) = pO + fer). (78) 


V 


Thus, if we define the fluid stress energy tensor as (71), with bọ = 1, we can see that the 
external forces only come from the dilaton field and the external electromagnetic field. This 
result is consistent with the one in [4]. It would be more clear to see the non-relativistic 
results by neglecting higher order terms of the solutions in the covariant form. The solutions 
n (63) become 


2 2 


ds? = + 2drdr + a l-V (r) dr? + 6,jda' da’ ] — 2u,dx'dr — 25 [1 — V(r)] vjda*dr 
2 
+ v'drdr + zB [rV'(r) Pdr? +(1-—V(r)) (v*dr? + vjvj;dz'dx') | 
r 2 2 2 igj (79) 
+ 7z [¢ K(r)Ogdr + (UF(r)oij + £ Fa(r)(o¢)is) dada | 7 
A=[|U(r)(1+ P+ v?/2) —rPU'(r)+ Ae) dr + [Uryu + Am) dz’, 
D =9(x) + LF(r) (O-g + vdih) + L Flr) (03H), 
up to €?, where the shear tensors are 
Oi; Oi; 
Tij = DGV) — a (0¢)ij = O@bOj)o — ras (80) 


and 0 = v, 0g = 670;¢0;¢ . Our aim is to obtain the forced NS equations, which turn 
out to be the constraint equations of the gravitational field equations (54) at order €. The 
solutions up to e? are enough to provide the forced NS equations, because higher order 
corrections do not make contribution in this order. With the solutions (79), the constraint 
equations at order €? give the incompressible condition wod;v' = 0, and at order e? give us 
with 

wo (rv; + 1A;v; + AP) — 8v = f + SO. (81) 


l 


These equations are just the temporal and spatial components of the equations in (75) up 
to e°. Here the external forces only have the spatial components as 


A) = -xo (07954318) + €o (403A), SP = -no (E? +BY), (82 


where B09 = F% and B” = a are the background electric and magnetic fields 
respectively. Further we define the dynamic viscosity and the normalized forces as 


Viw] = no/wWo, K = f Jo N. = FL? /wo, (83) 
the forced incompressible NS equations then become 
0,0; + 40,0; + AP — Up Po; = fO.4+ 6%, a=. (84) 


[w]i [w]i? 
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If we consider the characteristic scale L ~ e~! and the velocity v = yuv? ~ e, the Reynolds 
number of the dual fluid Re = vL/Vjy ~ 4aTy [1 + (nouo)/(Tuso)]. 

In addition, we can divide the electromagnetic forces into two parts fo = fre + fo, 
where 


flO) =— 9 (BO + IB), f= bonomo tamt OP), (85) 


where ge = Fe) and p = i Then p is just the Lorentz force from the extra 
electromagnetic field, while for from the boundary chemical potential. We now move the 
term i to the left-hand side of the NS equations. Defining the dynamic shear viscosity 
and the force density as 


m=z Doda fo - (86) 
HSo 4 Ty TH 80 Ty80 
we can rewrite the incompressible charged NS equations as 

3v + vavi + OP — vbu = FG) + AGP, dv! = 0. (87) 
In this case, the Reynolds number becomes Re = vL/vi] x Ty, proportional to the 


temperature of the fluid. 


4 Holographic forced fluid at cutoff surface 


In this section, we will generalize the previous discussions to the case of dual fluid at a 
finite cutoff surface by using the method which is introduced in Section 2. In this case, we 
only need to substitute 


gulr) = Ge (r) =r°V(r)/Ë,  Jaslr) = 1°/€°, (88) 


into the generic metric (1), and introduce a finite cutoff surface in the background solution 
and consider the non-relativistic expansions 


ð, ~ e, ĝ; ~ B ~ ðo ~ el, ð, ~ P ~ en (89) 
The equations of motion to be solved are given by (55) and (54). 


4.1 Charged fluid in non-relativistic limit 


Since we are considering the bulk solution with a finite cutoff in this section, following [23], 
we choose the gauge that gp, = 0, gra X Ŭa and g = 0. The perturbed solution up to 


ii 
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order e? with a Dirichlet boundary condition at the cutoff surface turns out to be 


ds? = — alee edly (%)dz°dF + s |- u — tha(Z) iy (Z) + Put) didi 
+ 7 KICA ae +r. F(r)/V (re) Fav + FA) Ea] dada", (90) 
A=A,@,2)dz" = Ess (7, #) + Ae (@)) B+ LA TE (2a) a) 
È =4(3°) + rË (r) VV re) (ëb) + 12 Fol(r) (ão) , 
where the gauge field relates to (63) through 
AGO (8) = -E il) Aa) = £4 Bi HAM ALE). (91) 


and n‘(Z) are given in (27) with a coordinate-dependent velocity. Again keep in mind that 
we have already taken the non-relativistic limit (89) of the solution (90) up to order e€? 
and higher order terms do not make contribution in the following calculations, they are 
reserved just for a covariant form. For the shear perturbations due to the boost, 


a a Th) ~ A x Pa A 
Foe f Say, Ga = Stn — 0 (92) 


and 6 = Ae Antty = P%9, ii,. For the shear perturbations due to the dilaton field, 


Š Po f" (yet — re~) - z 3 Pab 
For) = pei f Vy D (F4)ad = Aap — 7 


O5, (93) 


where 6. = P® (8.0050). The corresponding scalar perturbation equations give the fol- 
lowing solutions 


Ae g h(re) Jia g alre) E pP-l 
W)- Toes Wao): 


Te 
= l8 r?h(raV (r) ret 
K(r) = 2p(p = 1)r? i a r2 V(re) rp—1 


(94) 


Here two of the integration constants in (94) have been determined via the Dirichlet bound- 
ary condition of K(r) and Q(r). The other two integration constants a(re), hA(re), and the 
notation 

a(r) = 


e- o (95) 


V (re) r2 r?P 7 ret 
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will be determined via choosing the Landau gauge of the stress energy tensor of dual fluid 
and the induced charge current, whose expressions are given by 


Tas(Te) = (Kna — Kan) +E, FOr) = Nok. (96) 
The counter term at the cutoff surface is 
- 2 0 a ae 
(ct) x ap ~ 2 Nab 2 
755? = -tro [Erot r (aao - BG0")), (97) 


where ¢(r,) could be an arbitrary function of re. In fact, at the finite cutoff surface, the 
counter term is not necessary. To match the results at the AdS boundary in the case 
without the cutoff surface, we here add the counter term and can take č(re) = 1. 

The stress energy tensor and current in the derivative expansion of the fluid dual to 
the perturbed solution (90) can be written as 


a 


(98) 


Talr) = T+ Tg + TQ t+, TO = W(re)tiatin + Bre) thas, 
D (ra) = J® a IJO H -3 JO =F Ts) Ua 


where w(r.) = (re) + p(r-). The zeroth order parts in derivative expansion are just the 
stress energy tensor and current of ideal charged fluid. The energy density p(r-), pressure 
p(r-), and charge density (re) up to order e? in the non-relativistic expansion are given 
by 


P(Te) = Pa (re) — (reP)p, (ra); Po(re) = a Tr a 
LIGA 7. PolTe); (99) 


Alre) = Ng (Te) — (TeP)A,(Te), —-Mo(Te) = U"(re). 


P(e) = py (rc) = (r-P) py, (re) Drs) = 


The higher order terms in derivative expansion (98) are the dissipative parts of the 
p and 7 (g iss) Then using the 


stress energy tensor and current, we denote them by A 


Landau gauge up to e’, 
uF) = 0 > h(re) = Are), a F*9) = 0 = alr.) = hre), (100) 


so that g(r) = 0, and when re — oo we can recover the results (63) at the AdS boundary 
via transforming them to the coordinates in the boundary. - - 
After imposing the Landau gauge (100) up to order €?, we have GL) = JG =0 and 


TY = 28, (re)oa, T = —2ifg(r-)(Fs) av — Cp ('-)OsPav, (101) 


where the cutoff-dependent shear viscosity and entropy density are 


P p ed 
~ Th ~ 1 Th folre) 1 
e) =p = TFET: 102 
ir) = 3 Wigro > Fk (102) 
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And the coefficients associated with dilaton ¢ are 


T <_< Th alr) VJ- 
wd = ra [Bar + lr 1], 
x 4 Elfar V (Tre) (1.03) 
ro = —— |1 ar.) JV (re) — 
HO T E ve 
On the other hand, at the zeroth order, from (7) and (99), one has 
Pir) r?V' (rn) = tp.) — & bp hes) rv (re) 
Ar) = EE, tule) = Blo) tA) = EE, 00 


as well as ñ (rc) in (99), and S,(r-) in (102), thus one can show that the following thermo- 
dynamical relation still holds for the fluid at the cutoff surface 
A os 2 7 = 4 U(r.) = U(rp) 
Blr) = (re) -oA Alr) = LACIE a5) 
Te A/V (Te) 
where the chemical potential ñ, (re) is defined as the difference of the gauge potential be- 
tween the horizon and the cutoff surface. The dimensionless coordinate invariant diffusivity 
defined in [5] is 


5 fn) Pin felt) — 2 fi, Belraatre)] 
Dard = BoE E = ie Et Reged a 


which shows the dependence of the cutoff surface when the chemical potential is present. 


4.2 Forced Navier-Stokes equations 


As what has been done in section 3.1, we can redefine the stress energy tensor of the dual 
fluid as 
8 = W(Te)Uatp T Daa = 27 (Te) Fab, (107) 


by moving out the part of dilaton field. By using the constraint equations at the cutoff 
surface for gravitational field (72) and Maxwell field, up to e? we can get the forced NS 
equations and the conservation equation of charge current 


Pe =R Ot Gy PIO Har) Pt, =0, (108) 
where the external force from the dilaton field is 


r? 


A? (Te) = Ko lr.) (ioo) + Elre) (főpo), olr) = lr) = F, 


p= ~ 1 (109) 
pje E G La ate.) WG _ N ra 
ce p(p—1)/V(re) |r? (p— 2) (: re) VV (Pe) = V (re) i l 


21 


And the Lorentz force of the charged fluid is 
Prisi +h, RP =-ToFa RP =-ToFs, (110) 


with the external Maxwell tensor F es) = 23A” (re), and the induced Maxwell tensor 
pe = 20.40 (r-), where the gauge potential on the cutoff surface is 
LU (re) 


ACD (r) TA it (re) = (r-P) it, (re) | ta, Mo (re) = rey Vr) 


(111) 
If we furthermore redefine the following stress energy tensor by moving out the terms 
associated with the boundary chemical potential 


= = LG. ¢)So(Fe)thatl, + [By (Fe) — fo (Te)iba(Fe)| Mav — Zijo (Te) Sab- (112) 


where ře = r-(1 — P), when bo = 1, one has U(r) = 0, and ña (re) is just the chemical 
potential ñ, (re) of the fluid at the cutoff surface, the following relation holds up to order 


e, 


zaradu zfs) as F b Aa4-(s f plex 
F(T) — TE hig (Te) = Alr AP? MS ETE = f+ F. (113) 


This keeps the same form as in the case at AdS boundary discussed in the previous sec- 
tion. In this case, the external forces only come from the dilaton field and the external 
electromagnetic field. To see the above results more clearly, we write the them up to the 
desired order. For example, the solutions (90) up to order e° are, 


Vir) 
V (re) 


2 
d3 = + ae E + — 5|- V (r) 


Tre /V (re) V (re) 


pE A 
a7 - WD al B,dx'dr + (P+ 7 a £) drdr + H(r \( 28rd 
(re 


ave (veg ~via ) Pe ag Tay) (ae + ttt) 


(r265)drdr + 5 [rot )VV (re) +r Fol(r)(Go)is| da'de’, (114) 


, W(r) 7 — pav ia ~ _ Ulf) a zi 

E Ve) ral D (.- U(r ee ee 
Q(r) r? z 

a Sr E 

baserer r rore ) (5.o + pão) + r2F5(r)(5'5,6), 


and the shear components are given by 


dr? + ne = - h — | 2B,dx'dr 


7 a Ôi; X 7 x S nD 
Dij = dabi = = On (Oa )ag = dupo = > (115) 


where ĝ = 640,8;, 65 = 6%0,00;. At order e, the temporal component of the forced NS 
equations in (108) turns out to be 


w,(r.)0,8' = 0, w (re) = TeV" (re) /(E/V (re), (116) 


which leads to the incompressible condition, and the spatial component at order e? is 


tig (Te) (Ö-B: + Bibi) — reBy(reJAP — Mo (re)PBi = f(r) + Mee), (17) 
where only the spatial components of external forces remain 
FP (re) =- Kole) (878;95:0) + Elre) (595,06) , (118) 
FP (re) = — Fig(re) (FE + BEL) - (re) (FRO + RY). (119) 
If further redefine 
By = TPO) Bg = Mol). jo) = FO (re). jO- fr) aw) 
Wo (Te) Wo (Te) Wo (Te) Wo (Te) 
we can obtain the incompressible forced NS equations as 
0,8; + 81/8; + OP — Pð b = FS) + A, D=o. (121) 


Consider the characteristic scale of perturbations L ~ e~t and velocity 8B = ,/;G" ~ e€, 
the Reynolds number of the fluid turns out to be 
x BL 1 T (re) 


B= Fey Die) 


Lae fig (Te) Mo (Te) 
ae (re) h + F aol À (122) 


We can see that for uncharged black brane where the chemical potential vanishes, the 
Reynolds number of the dual fluid is proportional to the local temperature at the cutoff 
surface. Thus, when the cutoff surface approaches the horizon, the local temperature as 
well as the Reynolds number become larger and larger, and the fluid may become unstable. 
This instability may relate to the superluminal hydrodynamic sound modes when the cutoff 
surface is sufficiently close to the horizon [23, 44]. In addition, as in section 3.2, we can 


divide the electromagnetic forces into two parts: fo (r) = fom (re) + fev (re), where 


HO) = ilro) (AP -BAP tA — vA), 023) 
FPP (re) = folre) [fig (Pe) (3-8: + 618b) — rÑ, (re) OP] . (124) 
Define 
pf, = ee Tel FO FO (re) 7 = Fr) as) 
T (re)So (re) Ty (re)So (re) T, (re)So (re) 


and the dynamical shear viscosity Dj. = 7,(T)/[T, (1c) (Te)], the forced NS equations (121) 
then become 


8, Bi + 57058; + EPa — OO? B: = ie gi ee gp = 0. (126) 


In this case, the external forces only come from the additional electromagnetic field and 


dilaton field. And the Reynolds number has the form Re(re) œ Di (re) = AnT,(r,). 


5 Conclusions 


We have studied the thermodynamics and non-relativistic hydrodynamics of the holo- 
graphic fluid at a finite cutoff surface. As a calculation example, we have considered the 
case with the Gauss-Bonnet gravity. The isentropic flow and shear viscosity of the dual 
fluid have been obtained. The radial Einstein equation implies the isentropy of RG flow 
was first proposed in [5], and here we have generalized the discussion to the case of the 
Gauss-Bonnet gravity. The isentropy of RG flow can also be considered as an adiabatic 
process of the dual fluid. Note that instead of the entropy associated with the hologra- 
phy screen [45, 46], here the total entropy of the dual fluid is the horizon entropy of the 
background black brane, hence cutoff-independent. We have given a general formula (50) 
on the ratio of shear viscosity over entropy density of the fluid dual to the Gauss-Bonnet 
gravity. It shows that the ratio is independent of the cutoff surface. Namely it does not 
run with the cutoff surface. This may explain why the membrane paradigm [47, 48] and 
the AdS/CFT correspondence [49, 50, 51, 52, 53, 54, 55, 56] give the same result. This 
formula is also valid for the fluid dual to the Rindler spacetime in Gauss-Bonnet gravity. 

The main goal of this paper is to give the expressions of external force, coming from 
the bulk matters, of holographic fluid in the non-relativistic limit. To this end, we have 
considered an Einstein-Maxwell-dilaton system with a negative cosmological constant. By 
using the non-relativistic fluid expansion method, we have solved the system up to the 
second order of non-relativistic fluid expansion parameter €, and obtained the incompress- 
ible forced NS equations of the dual fluid at the AdS boundary and at a finite cutoff 
surface, respectively. The concrete expressions of external forces from the dilaton field and 
Maxwell field have been given. Here we have taken an new scaling of the dilaton field 
0,¢ ~ el in the non-relativistic limit, so that the external force provided by the dilaton 
field JP ~ e, meeting the scaling symmetry of the forced NS equations. Actually, in the 
derivative expansion of the stress energy tensor, these terms such as 0;60;@ are the second 
order dissipative terms of the dual fluid [3], and such terms may appear in the superfluid 
components [57]. However, in the non-relativistic limit, we move them to the right-hand 
side of the NS equations as external force terms. Note that here we have considered the 
case of minimal coupling of dilaton field. It would be interesting to extend this study to 
the case with non-minimal coupling [58, 59, 60). 

It turns out the Reynolds number of dual fluid is proportional to the local temperature 
of the cutoff surface in the uncharged case. Thus when the cutoff surface approaches to the 
event horizon of the black brane background, the local temperature and thus the Reynolds 
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number become larger and larger. This indicates that the dual fluid will become unstable 
when the cutoff surface is close enough to the event horizon. It would be interesting to 
further study the stationary turbulence by using the forced NS equations derived in this 
paper. In [61, 62, 20, 63, 64], the problems of turbulence have been studied via a holographic 
description with gravity. It was also shown that the nonlinear evolution of anti-de Sitter 
space might be unstable, and the energy of perturbations would be transformed to smaller 
and smaller scales like the turbulence energy cascades [65, 66]. Thus it would be interesting 
to establish the holographic turbulence using the AdS/CFT correspondence [67, 68]. 
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